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ABSTRACT 

We consider gauge theories with multitrace deformations in the context of certain 
AdS/CFT models with explicit breaking of conformal symmetry and supersymmetry. In 
particular, we study the standard four-dimensional confining model based on the D4- 
brane metric at finite temperature. We work in the self-consistent Hartree approximation, 
which becomes exact in the large- N limit and is equivalent to the AdS/CFT multitrace 
prescription that has been proposed in the literature. We show that generic multitrace 
perturbations have important effects on the phase structure of these models. Most notably 
they can induce new types of large- N first-order phase transitions. 
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1. Introduction 



In 't Hooft's large- N limit of gauge theories JIJ, the scaling of the bare gauge cou- 
pling g 2 ~ 1/N is tuned so that the vacuum energy is proportional to iV 2 . This scaling 
generalizes to an arbitrary action according to the rule: 

S = N 2 W(0 1 ,0 2 ,...) , (1.1) 

where W is a general functional of operators of the symbolic form 

O n = ^TvF n , (1.2) 

the set of single-trace gauge-invariant operators with expectation values of 0(1) in the 
large-iV limit!. For more general theories, including scalar fields and fermions in the 
adjoint representation, we extend the basic family of gauge-invariant operators to include 
these fields as well. These operators become quasi-classical in the large- N limit, in the 
sense that 

lim (0 0')= lim ( O ) ( O' ) . 

N^oo N— >oo 

This means that there is a notion of saddle-point configuration —a "master field" defined up 
to gauge transformations, which makes the 1/N expansion into a semiclassical expansion 
i- 

Known or conjectured master fields are usually established for single-trace actions, i.e. 
for linear W in (|1 . 1| ) , such as the Yang-Mills action. However, the behaviour of master 
fields under perturbations by multitrace operators is of primary interest, especially in 
the context of the AdS/CFT correspondence 0. In the holographic mapping, multitrace 
operators are associated to multiparticle states in the bulk theory. Hence they correspond 
to exotic deformations of the string background Moreover, truly non-perturbative 
effects in the bulk theory manifest themselves as finite- iV multitrace effects on the CFT. 
This is simply the translation of the fact that only O(N) elementary powers of the form 
Tr F n are algebraically independent: for n ^> N the single-trace operator decomposes as 
a sum of products of lower-order single-trace operators. Hence, the spectrum of the bulk 
theory must deviate significantly from a Fock space for states with O(N) "particles". 

It is then very interesting to study the effect of multitrace deformations on the 
AdS/CFT saddle point, particularly the effect of deformations that are non-polynomial 



2 We shall not discuss here operators with anomalous large- iV scaling. The most notable 
example is the theta-term with scaling (Tr F A F) = O(l). 
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in the traces. Recently, the AdS/CFT algorithm was modified to incorporate multitrace 
operators ]5],|6|] (see also J7|). Here we elaborate on some points made in || to argue 
that this modification can be understood in rather general terms, as an application of the 
mean-field approximation. 

In analysing the large- N master field, we could attempt a saddle-point approximation 
once we have managed to exactly integrate out 0(N 2 ) degrees of freedom. If we remain 
with O(N) degrees of freedom, this sets the order of magnitude of the fluctuations. Since 
the action is of 0(iV 2 ), we have a sharp saddle point. In practice, such a program only 
works in very restricted models in low dimensions, where we can integrate out explicitly 
the 0(N 2 ) angular variables (for a discussion of multitrace operators in these models, see 
||). Still, one can argue in great generality that in the leading large- N approximation W 
can be taken essentially linear. 

Let us suppose that we have managed to change variables in the path integral from the 
gauge field to the set of gauge-invariant monomials O n with n < O(N). In the process 
we generate a complicated (non-local) effective action T. At the large-iV saddle-point we 
have: 

dT dW 

<O cl ) + i^-(O cl ) = 0, (1.3) 



where we have incorporated the fact that the solution of the saddle-point equations is 
nothing but the planar expectation values: O c \ = limjv^oo^n)- 

In view of ( |1.3| ), it is clear that these equations are exactly the same as those that 
follow from a model with a single-trace action given by 

W(0) =J2(nO n , (1.4) 

71 

where the effective single-trace couplings Q n are given by 

dW 

C = m {OJ> • (1.5) 

Therefore, provided we only consider the planar N — > oo limit, any quantity of the orig- 
inal theory ( |1 . 1| ) can be computed in the single-trace theory ( |1.4| ), with the expectation 
values (O) being determined self-consistentlyJl Thus, in the AdS/CFT set-up, the com- 
bination dW(O c i)/dO n plays the role of the source for the single-trace operator O n , and 
this precisely determines the boundary conditions proposed in [H||. 



3 This argument assumes some explicit regularization, so that the path integral measure is 
denned independently of the details of the action. 
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Our discussion shows that the basic phenomenon is more general than the particular 
AdS/CFT set-up. Namely, it is a general consequence of the fact that the Hartree (or 
Thomas-Fermi) approximation becomes exact in the large- N limit (see for example 0). 
In this limit, the interactions between the gauge-invariant variables O n can be substituted 
by the interaction of each variable with a collective mean field that must be determined 
sefl-consistently. 

We should emphasize that these rules are only valid in the strict N = oo limit. 
The 1/N corrections will alter the master equation ( |1.5| ) since the Hartree approximation 
obtains corrections. Equivalently, the AdS/CFT boundary conditions of [|5||J] will receive 
1/N corrections, in addition to the usual loop corrections in the bulk of AdS. 



2. Master Field Dynamics 

To be more specific, let us suppose that the deformations by a certain single-trace 
operator O: 

5 S = N 2 ( [ d d xO , (2.1) 



are under control, in the sense that we are able to compute the planar one-point function 
(O)^ as a function of £ and the other couplings of the Lagrangian. Then we can compute 
any planar expectation value of the more general theory with perturbation 

Af 2 J d d x^ d F{iT d ° O) , (2.2) 

where F is general function, \x is a mass scale and do is the scaling dimension of the 
operator O in the single-trace model. We simply do our calculations in the single-trace 
theory (|2.1| ) with perturbation 



SS = N 2 C J d d xO , (2.3) 
where £ is given self-consistently by the solution of the "master equation" : 



G(C) = C - H a ~ ao F' [i~ do (€))■£ =0, (2.4) 

where the prime denotes differentiation. In principle, we can give £ a space-time de- 
pendence so that ( |2.4|) becomes a functional equation for an effective source. Such a 
generalization is appropriate to compute correlation functions in the multitrace-deformed 
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theory. However, for the purposes of this paper we are only interested in vacuum proper- 
ties of the master field, i.e. we consider only condensates and effective couplings that are 
translationally invariant in R d . 

The "master equation" ( |2.4| ) implies that multitrace deformations whose single-trace 
"elementary" operator has a vanishing one-point function are equivalent (in the large- N 
limit) to single-trace deformations, i.e. a constant shift 

5C = /- d ° F'(0) 

of the coupling dual to the single-trace operator O. Hence, in order to have specifically 
new phenomena associated to multitrace deformations we need non- vanishing condensates. 
This means that the auxiliary single-trace model with perturbation ( |2.3| ) must break con- 
formal invariance either explicitly or spontaneously. Since the one-point function depends 
on the particular state that we are considering, it is plain that the physical properties of 
multiple-trace deformations have a strong dependence on the full physics of condensates 
of the associated single-trace model. 

We may take F as a non-polynomial function of single-trace operators. However, we 
implicitly treat the non-linear terms as a perturbation since the scaling dimensions do are 
defined with respect to the single-trace theory. At any rate, it is interesting to evaluate 
( P-4|) when the function F becomes non-polynomial. 

Our main observation in this paper is the following. The function G( £ ) may have a 
complicated structure, being non-linear in both £ and the couplings of the bare Lagrangian 
W. In particular, if G( £ ) has various nodes, we have a set of solutions { C Q } f° r a given 
fixed value of the microscopic couplings in W. In this case we must select the master field 
that dominates the large- N dynamics among the various solutions Q a . 

By analogy with similar situations in large- N physics we characterize the dominating 
master field by requiring that the partition function be maximized: 

1 



lim — r log Z w = max 

N— >oo N 



i im . ^72 lo § Z (C a )w 



(2.5) 



TV^oo AT 2 

Large- N phase transitions induced by the multitrace couplings will arise when the dom- 
inating zero of G( C ) changes discontinuously as a function of the microscopic couplings 
in W . These phase transitions will be characterized by a "latent heat" release of 0(N 2 ). 
Typically, Z( £ ) will be a monotonic function of £, so that a change of branch in ( |2.5| ) 
will require that the cardinality of the solution set { ( a } changes as a function of W. 

Although the phenomena described so far are expected to be rather general, we will 
illustrate them in a specific example in the context of the AdS/CFT correspondence. 
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3. Multitraces in Deformed QCD 



As a concrete example along the previous lines, we consider a regularized version of 
four-dimensional non-supersymmetric Yang-Mills theory that has been introduced in fllCfl . 



In its most straightforward definition, the model is given by the low-energy theory on the 
world-volume of a stack of N parallel D4-branes at finite temperature T. Equivalently, we 
can view it as a Scherk-Schwarz compactification of the D4-branes on S 1 x R 4 , the compact 
circle having size 1/T. At large distances on R 4 the effective theory is a four-dimensional 
Yang-Mills theory modified at energies of 0(T) by remnants of the five-dimensional M = 4 
super Yang-Mills theory. The action is given by 



f <? x c = _L_ / d 4 x (ttf 2 + . . .) 

9ymJ 9ym J V 



(3.1) 



where the dots stand for other fields such as fermions and scalars of the parent M = 4 
theory, suppressed by powers of the cutoff scale \i = T. Planar quantities are functions 
of the 't Hooft coupling A, defined at the cutoff scale \x. In terms of the microscopic 
parameters of the parent D4-brane theory we have 

A = 5ym (aO n ~ 9s N fj, Va' , 

where g s is the string coupling and a' is the string's Regge slope. For A < 1 we have the 
standard planar perturbation theory of the four-dimensional Yang-Mills theory. On the 
other hand, for A > 1 we have a good description in terms of the low-curvature expansion 
of the black D4-brane metric. In this case, the expansion parameter is controlled by the 
curvature of the near-horizon metric in string units: 

a' 1 

m ~ a' 



with R c the curvature radius. Defining 



1 



the supergravity description is good for < x 1. At x ~ 1 we have the standard "cor- 



respondence point" in the sense of [JTJ , which represents the matching to the perturbative 
regime. As long as we only look at energy scales of 0(1) in the large-iV limit, we can 
neglect non-perturbative thresholds associated to large values of the dilaton, since these 
involve explicit powers of N. 
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The simplest multitrace perturbation in these models is a non-linear function of the 
Lagrangian density, 

S=^ J C + N 2 j fi 4 F(^- 4 C) . (3.2) 

According to ( |2.4| ), all physical quantities in this model, such as thermodynamic functions, 
condensates, Wilson loops, etc. , can be computed in the large- N limit in the auxiliary 
undeformed model 

- n 2 r 

S= T j£, (3.3) 
with effective 't Hooft coupling A given by the solution of the equation 

I" 1 = A- 1 +F / [a*-*</:> x ] , (3.4) 
where the gluon condensate (£)v is determined by the expectation value of the action: 

^)_=Vol(R*) (^^ + ...)_=i,A 2 |logZ(A). (3.5) 
The partition function in the planar supergravity approximation is defined in terms of the 



thermal free energy of the D4-brane (see, for example ||12|| ) 

1 



Vol(R 4 



log Z(X) = N C Xfi , (3.6) 



where C is a positive numerical constant. This expression for the partition function has 
been normalized to the Euclidean action of the wrapped D4-brane metric with supersym- 
metric boundary conditions, i.e. we define the five-dimensional thermal free energies with 
respect to the T = vacuum. 

Notice that, even if the general multitrace deformation of the T = D4-brane theory 
may break supersymmetry, the N = oo effective theory ( |3.3| ) does not. Hence, the D4- 
brane theory reduced on a supersymmetric circle will be supersymmetric at N = oo and 
no condensates will be induced .@ This implies that the condensates are entirely due to 
thermal effects of the D4-brane theory and our normalization of (|3.6|) is the physically 
correct one. 

Combining (|3.5|) and (|3.6| ) we find the value of the gluon condensate (c. f. |]T~3H ) : 

(C)j = CX 2 . (3.7) 



4 Casimir energies are not induced either, since the D4 world-volume is fiat S 1 x R 4 . 
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This expectation value has the crucial property of diverging as A — > oo. Since this is pre- 
cisely the supergravity regime of the effective single-trace theory, we learn that multitrace 
deformations are potentially stronger in the region where AdS/CFT is under quantitative 
control and they may be reliably studied. 

In terms of the dimensionless expansion parameters x = 1/A and x = 1/ A the master 
equation reads 

G{x) = x — x — F f [C/x 2 ] = 0. (3.8) 

Equation ( |3.8| ) was derived within the supergravity approximation to the near-horizon 
black D4-brane solution. In terms of the supergravity expansion parameter x, this is the 
regime: 

< x < 1 . (3.9) 

As before, these limits ignore other thresholds that are related to large dilaton corrections 
and are of subleading order in the 1/N expansion. 

One important property of (|3.8|) is the redundancy of the description in terms of the 
original variables in the microscopic Lagrangian, i.e. the coupling x and the multitrace 
couplings that define the function F. For a fixed value of x all models in the codimension-1 
submanifold 

: x-x + F' = 

have the same large-iV properties. The region of the microscopic coupling space where 
supergravity is a good approximation is the union of these submanifolds for < x < 1: 

S = |J M- x . (3.10) 

0<x<l 

One component of the boundary is Mx=o defined as 

M : x + F / (oo) = 0. (3.11) 

It yields the strong-coupling (low-curvature) limit of the AdS/CFT background. On the 
other hand, the correspondence line (the matching to perturbative variables) occurs at 
x = 1 or 

Mi : x + F'(C) -1 = 0. (3.12) 
Although MqUMi are components of the boundary of 5, they do not exhaust it in general. 
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3.1. Multicritical Behaviour 



For < x < 1 there is always a standard solution of ( |3.8j ) that is valid for very small 
multitrace couplings. This solution has x ~ x and can be obtained iteratively as the limit 
of the set {#(&)} with 

X( k +i) = x + F' [C/(x( fc )) 2 ] , X(q)=x. (3.13) 

However, it is clear that there will be other solutions if F'[C/x 2 ] shows "bumps" in the 
supergravity interval < x < 1. 

Let us assume that F' admits a finite Laurent expansion around the origin, so that 
the master equation takes the form: 

G(x; x,f J )=x-x-J2^k = - ( 3 - 14 ) 

The j = term is equivalent to a constant shift of x and has been removed from ( |3.14| ) . 

Our first result is a simple consequence of the divergence of (|3.7| ). The pole part of 
F, corresponding to j < in ( 3.14| ), has no dramatic effects in the supergravity interval 



< i < 1. Thus, multitrace deformations that are completely singular in perturbation 
theory become rather tame in the supergravity approximation. This looks surprising at 
first sight, but it fits naturally with the character of AdS/CFT as a strong/weak coupling 
duality with respect to the 't Hooft coupling. 

Conversely, perturbations that are polynomial in multitraces translate into non- 
analytic contributions to G(x) and therefore dominate the supergravity regime at x — > 0. 
In this limit G(x) diverges with a sign that is correlated with that of fj, J being the largest 
value of the index j. In particular, for fj<0 and small there is always a solution: 

This solution disappears for fj > 0, unless one also dials the microscopic 't Hooft coupling 
to negative values: x < 0. 

We have found that for < x < 1 and small \fj\, we have a discrete jump in the 
number of solutions of the master equation as fj crosses zero. This is a source of possible 
phase transitions. 
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A more general multicritical behaviour in the vicinity of x ~ will depend on the 
higher multitrace powers. Let us consider a simple example of a deformation proportional 
to 

| / ^ (TrF*) 2 + (Tr F 2 ) 3 . (3.16) 



The master equation (|3.14|) reads 



fi h 

G(x; x, fi, / 2 ) = x - x - ^2 - 34 = , 

with fi ~ gi up to numerical constants. Besides the standard solution x + ~ x for very- 
small fj there are other interesting solutions. Consider x > 0, fi > and /1 < 0, with 
$2 <C |/i I C 1 and furthermore <C /i- Then, the master equation has two small 
solutions in the vicinity of 

x- ~ yj-hlh, x '- ~ V~fi/ X ■ 
These solutions coincide for / 2 ~ x fi and disappear for larger values of fi- 

3.2. Phase Transitions 

The previous discontinuities in the solution set of the master equation translate into 
large- N phase transitions. Since the partition function scales as 

N 2 

logZ(x)oc — , (3.17) 

x 

we find that the dominant solutions in the supergravity approximation are those with the 
smallest value of x within the unit interval. The jump of the effective action across the 
transition from x a to xp is given by 



log 



Vol(R 



Z[Xg) 



N^Cfi 4 { — -—) ■ (3.18) 



Coming back to the examples in the previous subsection, we see that there is always a 
phase transition when fj crosses zero from negative to positive values. In this case x a = 
and x/3 ~ x > 0. The density of "latent heat" in (|3.18j ) is infinite. This phase transition 
is not hard to interpret. Since fj is the coupling of the multitrace interaction of highest 
order, it dominates the limit of large field-strengths. Hence, the very strong singularity for 
fj — > 0~ reflects the fact that the microscopic action is not bounded below for fj < 0. 
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A more physical phase transition with finite "latent heat" takes place in the two- 
coupling model (|3.16| ) with x > 0, > and /i < 0, when the two solutions around 
X- ~ y—fijx coalesce as we decrease the magnitude of I/1I//2. For small values of this 
ratio the only solution is x ~ x. 

This example illustrates the general pattern of phase transitions in this class of models. 
When the minimal solution x_ of the master equation is separated from the first subleading 
one x'_ by a local maximum of G{x) 1 a variation of the parameters can bring the maximum 
to zero and make the two solutions coalesce X- = x'_. A further variation of the parameters 
can bring the maximum to negative values and make the double solution disappear. This 
generic situation is depicted in Fig. 1 below. 



G 




Fig. 1: A depiction of a typical phase transition. The solid line shows the 
function G(x) with three zeros, x_ < x'_ < x + . The dotted line shows the 
degeneration of the lower zeros x_ = x'_ and their disappearance in favour of 
x+. When the partition function is dominated by the smallest solution this 
degeneration yields a large-N phase transition. 
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4. Generalization to Other Dimensions 



This set-up can be generalized to the regularized Yang-Mills model on R p , with p < 5, 
in terms of a hot Dp-brane model and the corresponding generalization of the AdS/CFT 
correspondence [0. In this case, the effective dimensionless 't Hooft coupling normalized 
at the cutoff scale [i = T is given by 

This is the expansion parameter of the planar perturbative expansion. The expansion 
parameter of the supergravity approximation that arises at X p (p) 3> 1 is: 

The large- N solution of these models perturbed by multitrace interactions of the form 

*j**„F^C) (4.2) 

can be studied along lines similar to the p = 4 case above. Here, C denotes the Yang-Mills 
Lagrangian operator, corrected by regularization artefacts at the scale \i = T '. As before, 
one reduces the problem to the study of an effective single-trace model with supergravity 
expansion parameter 

x=(l/X p (fi))^ (4.3) 

that is determined self-consistently. The supergravity regime of the N = oo problem is 
then given by < x <C 1. The partition function in the single-trace model with effective 
coupling x is 

' log Z( x ) = N 2 (5 - p) C p fi p x s ~ p , (4.4) 



Vol(R 

where C p is a positive numerical constant. The gluon condensate is given by 

^- 4 <0, = (p-3)%. (4.5) 

These expressions show that the p = 3 case, based on the hot D3-brane, yields trivial 
multitrace deformations in this approximation. This is a consequence of the free energy of 
D3-branes being very smooth for large 't Hooft coupling. Of course, this situation changes 
when considering subleading terms in the a' expansion of the supergravity background. It 
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is interesting to study these corrections in more detail, although we will not attempt to do 
this here. 

For p < 3 one finds a situation somewhat similar to that discussed before in the p = 4 
case. The master equation for x reads: 

x 5 -P = x 5 - p + F'[(p-3)C p /x 2 ] = 0. (4.6) 

Hence, the same qualitative properties follow, regarding the multiplicity of solutions at 



small x. In particular, the crucial singularity at x = of the gluon condensate (|4.5|) still 
holds. 



The main difference with p = 4 is that, according to (|4.4|) , for p < 3 it is the largest 
solution x+ that dominates the partition function. As a result, we expect that the standard 
solution x ~ x will dominate and that sharp phase transitions will be more difficult to 
produce than in the p = 4 case. 



5. Conclusions 



In this paper we have studied some simple multitrace deformations of the basic non- 
supersymmetric QCD model in [ID], as well as its generalizations to less than four di- 
mensions. In particular we have considered deformations by a non-linear function of the 
Lagrangian operator. 

Our main result is the emergence of new types of "multicritical" behaviour, similar 
in many ways to those studied in the context of matrix models ||. There appear various 
competing master fields whose dynamics yields new examples of large- N phase transitions. 
It turns out that the dynamical effect of multitrace deformations is particularly strong in 
the supergravity approximation to the AdS/CFT master field. 

These results suggest various avenues for further research. It would be interesting 
to study more examples of large- N phase transitions induced by multitrace deformations. 
Eventually, these phase transitions should be related to the breakdown of string perturba- 
tion theory in the geometrical description of the large- N master field. Another interesting 
question is the effect of multitrace deformations on other large- N phase transitions that 
have been identified in single-trace models, in particular, the phase transitions associated 
to theta-dependence in [I5[] or those related to finite-size effects, as in [|ll^|10| , |T2] . 
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